QUANTUM ERGODICITY FOR RESTRICTIONS TO 

HYPERSURFACES 



SEMYON DYATLOV AND MACIEJ ZWORSKI 

Abstract. Quantum ergodicity theorem states that for quantum systems with er- 
godic classical flows, eigenstates are, in average, uniformly distributed on energy 
surfaces. We show that if is a hypersurface in the position space satisfying a 
simple dynamical condition, the restrictions of eigenstates to A^ are also quantum 
ergo die. 



1. Introduction 

In a recent paper [13] Toth and Zelditch proved a remarkable result stating that if 
(M, g) is a compact manifold with an ergodic geodesic flow, then quantum ergodicity 
holds for restrictions of eigenfunctions to hypersurfaces satisfying a certain dynamical 
condition. In an earlier paper [12] they established a similar result for bounded domains 
in whose boundaries are piecewise smooth and whose billiard flows are ergodic. 

The purpose of this note is to provide a short proof of a semiclassical theorem which 
simultaneously generalizes both results. Our approach avoids global constructions 
and calculations by reducing equidistribution for restrictions to equidistribution in 
the ambient manifold. The geometric condition (1.7) enters to obtain a decorrelation 
between contributions to the restrictions coming from different parts of phase space. 
The proof uses some ideas of [5, Appendix D] but we do not refer to any results from 
that paper. 

For the standard quantum ergodicity result established by Shnirelman, Zelditch 
and Colin de Verdiere, see [5], [7], [12], [13] and references given there. The case of 
Riemannian manifolds with piecewise smooth boundaries was established in a special 
case by Gerard-Leichtnam [6], and by Zelditch and the second author [14] in general. 
A semiclassical version of quantum ergodicity was first provided by Helffer-Martinez- 
Robert [7]. 

To make the presentation more clear, in the introduction we will work in the setting 
of manifolds without boundary, referring to Appendix A for modifications in the case 
of manifolds with piecewise smooth boundaries. 

Let (M, g) be a compact smooth Riemannian manifold and consider 

P{h) := -h^Ag + V{x), VeC°^{M;R). (1.1) 
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Figure 1. Left: the situation prohibited almost everywhere by the dy- 
namical assumption (1.7). Right: The projection map tte and the re- 
flection map jE in the cotangent space over some point x & N. 

The operator P{h) is self-adjoint when acting on half-densities (see [15, Chapter 9]), 
L'^{M,Q\^j^). (This technical point is helpful when more general operators are consid- 
ered.) The classical symbol of P{h) is given by 

p{x,0 = \^\l + Vix), ix,OeT*M, 

and p defines the Hamiltonian flow, 

:= exp{tHp) : p-\E) p~\E), E eR, (1.2) 

We make the following assumption on a range on energies: 

For E G [a, b], dp\p-i(^E) 0; ^^d the flow ipt : p^^{E) — )• p^^{E) is ergodic, (1.3) 

where ergodicity is with respect to the Liouville measure /i^ on p^^{E). 
Now, let iV be a smooth open hypersurface in M. We define 

EE:=p'\E)n7T-\N), (1.4) 

where vr : T*M — )■ M is the natural projection. We note that S^; is a smooth hyper- 
surface in p~^{E) if 

V{x) = E ^ dV{x) i N*N, (1.5) 

and for simplicity we make this assumption for all E G [a, b]. For > it is satisfied 
when V = 0, and that is the setting of Theorem 3. 

By restricting elements of S^; to TN we obtain a map 

TiE-.^E-^ Be := vrs(Ss) C T*N, (1.6) 

which is a local diffeomorphism almost everywhere. It defines a unique nontrivial 
involution 
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which is the reflection across the orthogonal complement of the normal bundle A^*A^ C 
T^M with respect to the metric on the fibers of T*M induced by g. This involution 
enters into the dynamical assumption we make on (see Fig. 1): 

For E G [a, b], the set of p G S^; satisfying ipt{p) G S^; 

(1.7) 

and ftilEip)) = lE{ft{p)) for some t 7^ 0, has measure 0. 
We denote by Uj{h) a normalized eigenfunction of P{h) with an eigenvalue Ej{h), 
{P{h)-E,{h))u,{h)=0, ||«,(/.)|| . =1. 

To formulate the quantum ergodicity theorem for restrictions, we need to restrict 
half-densities to and that requires making a choice. Suppose / G C°°{M), /Itv = 0, 
df\N 7^ 0. Informally, the restriction is now defined using, \dx\2 = \dy\2\df\^ , x G M, 
y E N. More precisely if, in local coordinates, x = (x',x„), = {xn = 0} then, in the 
half-density notation of [15, §9.1], 



(u{x)\dx\2) Itv := u{x', 0)\dx'\^ 



Using the notation of [15, Chapter 14], reviewed in Section 2 below, we can state 
our main result. See Appendix A for the modifications needed in the case when M 
has a boundary. 

Theorem 1. Suppose that {M,g) is a compact Riemannian manifold with a piecewise 
smooth boundary satisfying (A.l) and that uj = Uj{h) are normalized eigenfuctions of 
the Dirichlet realization of P{h). 

Suppose also that (1-3) holds and that N is a smooth open hypersurface not inter- 

1 /2 

secting dM and satisfying (1.7). For Q G \E'™(M, ) put vj := Quj{h)\N, where the 
restriction operator on half densities is defined in (1.8). Then for A G \&°(A^, fi],^^), 
compactly supported in N, we have 

^" 5Z (^^^■'^i)L2(iv,<^)- / T^%cT{A)\a{Q)\' duE, 

where (j{A) G S'^{T*N) is the symbol of A, a{Q) G S"'{T*M) is the symbol ofQ, and 
the measure ve is defined onY^E by the identity 

—^%^,=duEAdf, (1.10) 
with pe the Liouville measure and f defining the restriction of half- densities in (1.8). 

Remarks, (i) The measure ue has a particularly nice description in the case V = 0. 
Assume that E = 1, then Si = S%M where S^jM C S*M denotes the cosphere bundle 



0, h^O, (1.9) 
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of M restricted to A^. The Liouville measure, /ii, on S*M, induces, for each x G M, a 
measure on S*M, fj,^, such that 

= [ fx^inn s*M)dvo\g{x), n c s*m. 

Our measure on Si = S*^M is then given by 



''i(^) = (ci*M\ / -"-(^ ^ 5:M) rfvol,|,(a;), T C S^M, (1.11) 

where g\N is the metric on N induced by g. (Here Vt and F are Borel sets.) See 
Appendix B for details. 

(ii) The now standard argument due to Cohn de Verdiere and Zelditch and described 
in [15, Theorem 15.5] shows that this resuh provides pointwise convergence for a 
density one subsequence: there exists a family of sets, A(/i) C {a < Ej < 6}, such that 
lim/i_>o ^{f<Es<b} ~ 1' ^'^•^5 notation of Theorem 1, 



sup 

EjeA{h) 



/ 'n%a{A)\am' duE, 



0, as/i^O. (1.12) 



In step 4 of the proof of [15, Theorem 15.5], to pass from (1.12) for a countable dense 
family Ak of pseudodifferential operators to the full statement, one needs uniform 
boundedness of Vj on L'^{N). While this need not be true for the whole sequence, it 
holds for Ej G A(/i) if we take one of the operators A/^ to be the identity. 

(iii) The dynamical condition of Toth-Zelditch [13] is stated using Poincare return 
times but the analysis in that paper shows that it is equivalent to our condition (1.7). 
The paper [13] provides interesting examples for which it is satisfied. 

(iv) We sum over eigenvalues in a fixed size interval [a, b] in (1.9) since the correspond- 
ing smoothed out spectral projectors are pseudodifferential operators. It would be 
interesting to prove an analogous statement for size h intervals, whose projectors have 
more complicated microlocal structure, as in [5, Appendix D]. 

The paper is organized as follows. In Section 2 we review basic concepts of semiclas- 
sical quantization and present slightly non-standard facts needed in the proof. The key 
point is that even in the case of manifolds with boundary we only need to work with 
standard pseudodifferential operators. In Section 3 we present a general decorrelation 
result and in Section 4 the proof of Theorem 1. Since no reference for quantum er- 
godicity for semiclassical boundary value problems seems to be available. Appendix A 
present a proof in the spirit of [14] with some simplifications based on [5, Appendix D] 
(see Lemma 2.2). The high energy result for Laplacians, as presented in [13], follows 
from Theorem 1 but there is an issue at energy 0, and Appendix B shows how that is 
overcome. 
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Except for the efficiency provided by direct semiclassical methods [15], the proofs 
are similar to those in [13]. The one significant difference is the treatment of quantum 
ergodicity for microlocal Cauchy data in Section 4 - see [3] for comparison. 

Acknowledgements. We would like to thank Nicolas Burq, Oran Gannot and 
Stephane Nonnenmacher for helpful comments on the ffist version of this paper, Hart 
Smith for explaining an alternative proof of Lemma 2.1 based on finite speed of propa- 
gation, Steve Zelditch for encouraging us to handle the boundary case, and an anony- 
mous referee for many helpful suggestions. The partial support by National Science 
Foundation under the grant DMS-1201417 is also gratefully acknowledged. 

2. Semiclassical preliminaries 

We will use the calculus of semiclassical pseudodifferential operators described in [15, 
§9.3, §14.2]. Our operators will always be supported away from the boundary of the 
manifold and hence can be considered as operators on a boundaryless manifold (for 
example, by considering the double space of a neighborhood of our manifold in a slightly 
bigger open manifold without boundary). A notable exception is the Schrodinger 
operator —h'^Ag + V{x), which can however be extended smoothly past the boundary. 

For a compact manifold, X (which could be different from the compact manifold M 
considered in Section 1), the class ^'^(A) denotes operators of order m, so that, for 
instance —h'^Ag G \E'^(M). We have the symbol map, cr, appearing in the following 
exact sequence 

/i^;^-^(A) — ^ ^^(A)^ S'^{T*X)/hS'^-\T*X) — ^ 0, 

where S"^ denotes the standard space of symbols. We take some quantization map 
Op^ : 5™(T*A) ^ ^]!^; it satisfies 

a{Opf,{a)) = a mod hS"'~\T*X). 

We also introduce the class of compactly microlocalized pseudodifferential operators, 
^™'"P(A): A e ^-°°(A) is in ^^^(A) if for some x e C,°°(T*A), 

Op,ii-x)Aeh^^-^{X). 

For this class the definition of WFh{A) given in [15, §8.4] applies. From the same 
section we take the definition of microlocal equality of operators. 

Following [4, §2.3], [10, §3], and [15, §11.2] we consider Fourier integral operators 
quantizing a canonical transformation k : f/i — )■ U2, Ui (s T*A and U2 (s T*Y, k 
defined on a neighbourhood of Ui: we say that an operator F : L'^{X) — > L'^{Y), 
quantizes k if for any A G ^il™{Y) with WFh{A) d U2, 

F*AF = B, Be ^pP(A), a{B) = K*a{A). (2.1) 
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We further require that F be microlocally unitary in the sense that = F* mi- 
crolocally near Ui x f/2. If F quantizes k, then the operator F* quantizes k~^. More 
generally we say that G is a Fourier integral operator associated to the canonical 
relation k, if G = AF, for F above and some A G \1/°(X). 

The standard example is given by F(t) = Q-^tPW/^^ where P{h) = —h'^Ag + V{x) E 
^^(M) (or a more general operator) which quantizes the Hamiltonian flow ipt : = 
exp{tHp). 

We say that a tempered operator (see [15, §8.4]) G : L^{X) L^(F), is compactly 
microlocalized if for some A G \l/™™^(y) and B G \E'™™''(X), 

AGB -Ge (2.2) 

In that case we can deflne WF'^(G) C T*Y x T*X, by taking the twisted WFh of its 
Schwartz kernel, Kq: 

WFUG) := {{y, rj; x, : {v, ^, e WF^iKc)}. (2.3) 

If F is associated to some canonical transformation k, then WF'f^{F) lies inside the 
graph of K. 

We recall from [15, Theorem 14.9] that if M has no boundary, and P{h) = —h'^A + 
V{x), f e G^{R), then 

f{P{h)) G v]/— P(M), a{f{P{h))) = f{p). (2.4) 

We now present three lemmas which will be used in the paper. We assume that the 
manifold M and operator P{h) are as in Appendix A. 

The flrst lemma makes an observation that f{P) is a nice operator away from the 
boundary. 

Lemma 2.1. Suppose that f G C~(M), x e C~(M°). Then 

f{P{h))x = Af + OmM)^L^M){h^), (2.5) 

where Af G \1/^°™''(M°) is compactly supported away from the boundary and its principal 
symbol is given by f{p{x,C,))x{x). 

Proof. Let ipo G G^{{—6, 6)) be equal to 1 near 0. We write f{P)x using the Schrodinger 
propagator, 

fiP)x=7^ I f{t/h)e''''/\dt=-^ [ f{t/h)Mt)e''''/'xdt + OL^M)^LHM){hn, 

where the error estimate follows from the decay of / and the untarity of e**'^/'^. We 
then write 

f(P)x=^[ [ f{r)e'i^''-^hMt)dtdT. 
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Choose 6 sufficiently small so that for \t\ < 6, p{x,^) G supp /, and x G suppx, we 
have dg{7T{ipt{x,^)),dM) > e > 0. In that case we can use the local parametrix for 
gitp/h^ see for instance [15, §10.2]. An application of the stationary phase method in 
{t, t) variables gives the conclusion of the lemma. □ 

As before, let be the full orthonormal system of eigenfunctions of P{h) 

with eigenvalues Ej{h). Lemma 2.1 applied to the operator f{P{h))A, where A G 
^'™(M) is supported away from the boundary, and / G C^(]R), gives, together with 
the trace formula [15, Theorem 14.10] 

i2nhy^y2fiE,){Au,,u,)= [ fip)aiA)dfi^ + Oih), (2.6) 

J Jt'-m 

where fia is the symplectic measure, fi^ = a"'/n\. 

The second lemma, in the spirit of [5, Appendix D], gives estimates using norms 
of symbols: 

Lemma 2.2. Suppose that A G "^^{M") is a pseudodijferential operator compactly 
supported away from DM . Then for each a' < a < b < b' , 

{27chr Il^^.lli. < lk(^)lli.(p-i([a',6'])) + ^W (2-7) 

Ej€[a,b] 

where the norm of a {A) is taken with respect to the measure fi„. 

More generally, if N G M is a fixed smooth submanifold ( of any dimension ) not in- 
tersecting DM , then there exists a constant C such that for each A G \E'^(A^) supported 
in a fixed compact subset of N , 

J2 P(^.U)lli^<C'lk(I)|li.w,-i(K6q)nT*M)) + 0W (2.8) 

Eje[a,b] 

Here T^M is the cotangent bundle of M restricted to N and vr : T^M — )■ T*N is the 
projection. 

Remark. We note that in the case when A = 1 we recover the bound 

/i" J2 <C. (2.9) 

Eje[a,b] 

By contrast, for individual eigenfuctions the bound Ch^^ is optimal - see [1] and [11]. 
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Proof. To show (2.7), take / G C^{a',b') such that < / < 1 everywhere and / = 1 
on [a,b]. Then we write by (2.6), 

i27Thr Yl p«.iii^<(2^/^)"5^/(i5i)(A*^w„w,) 

= [ f{p)\a{A)\'dfi, + 0{h)< [ \a{A)\^dfi^ + 0{h). 

Jt*M Jp-^([a',b']) 

To show (2.8), denote by Rn '■ C°°{M) — > C°°{N) the restriction operator and note 
that 

Eje[a,b] j j 

= h-\\AR^m\\i,. 

The Hilbert-Schmidt norm on the right-hand side is equal to the norm of the 
Schwartz kernel K of AR]\ff{P). Recall that f{P) is a pseudodifferential operator 
when localized away from the boundary, by Lemma 2.1. Note that K = 0{h°°) away 
from the diagonal of embedded in x M. To estimate K near the diagonal, we 
choose local coordinates x = {x',x"), x' G M'^, x" G M""*^, where k = dimA^, on M 
near some point of N, in which N is given by {x" = 0}. If a is the full symbol of 
A in these coordinates (in the standard quantization) and b is the full symbol of the 
pseudodifferential operator f{P{h)), then we can write 

ir(x,?/) = (27r/i)-"-'= / eT^^''-'^-'-'^^'-^'-y^^d{x,ri)b{z,0,^)dzdrid^, 



here y,^ G M" and x,z,t] G M'^. By the unitarity of the (semiclassical) Fourier trans- 
form, the Ll y norm of K{x,y) is equal to the L^.^^ norm of 

iri(x,0 = (27r/i)-"/'-'^ / es("-''-^-''+^-«')a(x,r/)6(z,0,Orf^t^^- 



The method of stationary phase shows that 

Now, /i"/^ times the norm of Ki is bounded by a constant times the norm of a 
on the set 7r(supp b fl T^M), with an 0{h) remainder. □ 

To formulate the next lemma we define 

Diag(T*M) := {(p,p) : p G T*M} C T*M x T*M. 

Lemma 2.3. Suppose thatG : L'^{M) -> L'^{M) is a compactly microlocalized tempered 
operator in the sense of (2.2), compactly supported away from the boundary, and that 
f G (M) . Then for G satisfying 

WF'JG) n Diag(T*M) = 0, 
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we have 



(2.10) 



j 



Proof. The left-hand side of (2.10) is equal to the trace of Gf{P{h)). We can write 
G as a finite sum of operators of the form X1GX2, where Xi,X2 G satisfy 
WFfe(Xi) n WF;,(X2) = 0. Then by the cychcity of the trace, 



In the proof of Theorem 1 we will encounter expressions involving {Fuj,Uj), where 
{uj{h))jQ^ is the full orthonormal system of eigenfunctions of P{h) = —h'^Ag + V{x) 
with eigenvalues Ej{h), and F is a compactly microlocalized semiclassical Fourier 
integral operator. This section shows that the sum of such terms over j in an 0{1) sized 
spectral window is negligible when the canonical relation of F satisfies a 'nonreturning' 
assumption; we call this phenomenon decorrelation for Fourier integral operators. 

Assume F is a compactly microlocalized tempered operator L'^{M) -4 L'^{M), in 
the sense of (2.2), and 



where the wavefront set is defined in (2.3), k : — )■ V2 is a canonical transformation, 
Vi,V2 G T*M are open sets, and Kj C Vj are compact sets such that k{Ki) = (In 
our case, F will be a Fourier integral operator, but this is not required in the proof.) 

For each t G M, define the t-exceptional set. 



£^it):={peK^nip.t{Ki):ipMp)) = K{ift{p))}, ipt:=exp{tHp). (3.2) 



Here (ft is the broken Hamiltonian flow of p, defined almost everywhere and described 
in Appendix A. 

The decorrelation result is given as follows: 
Lemma 3.1. Suppose that a < b are fixed, and that there exists to > such that 



TT{XiGX2f{P)) = Tr(X2/(P)XiG) = 0(/i°°), 
as X2f{P)Xi e by Lemma 2.1. 



□ 



3. Decorrelation for Fourier integral operators 



FlU^^x. = 0(1), WF;(F) C {{k{p),p) : p G K,}, 



(3.1) 




(3.3) 



\t\>to 



where S^it) is given by (3.2) and p^ is the symplectic measure. 

Then for each F compactly supported inside M° and satisfying (3.1) 




(3.4) 



Ej<^[a,h] 
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Proof. Take T > to and denote 

KT:={{BTUK-\BT))nK,)u [j £^{t). (3.5) 

to<\t\<T 

Here Bt is the set of points at which the broken Hamiltonian flow is not well-defined 
on the interval [— T, T], see (A. 2). Then is a compact subset of Ui and ^^{Kx fl 
p~^{[a, b])) = 0. Therefore, there exists an open set Ut C Ui and constants a' < a and 
h' > h such that 

Kt C Ut, ^l,{UTr\p-\[a\h'])) < T'K 

Take Xt G ^™'^p(M) satisfying \a{XT)\ < 1, WF/,(Xt) C Ut, and Xt = 1 mi- 
crolocally near Kt- Since F is bounded on L'^{M), \{FXTUj,Uj)\ < CWXtUjUl^. 
Hence (2.7) and (A. 3) give 

J2 \{FXTU,,u,)\<c(h- Yl II^T^.-lli.)' 

Ej(^[a,b] ^ Ej&[a,b] ' (3-6) 

< C(||a(X^)||i.(^_.(f,, + 0^(/^))i < C(T-i + 0^(/.))i, 

where C denotes a constant independent of T and h. 

We now analyse the contribution of F\ := F{1 — Xt)- For that define 

^ Jo 

For each eigenfunction Uj, we have 

(Fl Uj ,Uj) = {{Fl) TUj ,Uj). 

We now take some / G C^(]R) such that < / < 1 everywhere and / = 1 near [a, b]- 
Then by (A.3), 



Ejia[a,b] Ej£[a,b] 



<C(/."5^/(F,)||(i^l)TU,|| 

We now write 



(3.7) 



j 3 

Jo Jo 
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Since \\e'^^'^~'^^^^'^^^'^ F*e'^^^~'^^^^'^^^^ Fi\\l2_:^i2 is bounded uniformly in t,s,h, we estimate 
the integral over the region |t — s| < to using the upper bound on the number of 
eigenvalues, (A.3), 

^ L, <T ^" 5Z /(^i)(e'^^-*^''^^^/'^Fi*e'(*-^)^('^)/'^FiM,-, uj) dtds < CT-\ (3.8) 

|t-J|<to j 

where C is again a constant independent of T and h. 

It remains to estimate the integral over the region to ^ l^"*! < T. For Iq < \r\ < T, 
define the operator Gr = e^'^^^'^^^^F^e~^'^^^^^^'^Fi, and rewrite the studied integral as 



h / h''y^ f{Ej){G,^tUj,Uj)dtds. 

. J 0<t,s<T ^—^ 



2^2 

" |t-s|>to 3 

Take x ^ C^(M°) such that F\ = Fix, and thus Gr = GrX, then by the cyclicity of 
the trace as in the proof of Lemma 2.3, we can replace Gg-t by xGs-t with an Orih'^) 
penalty. However, by Lemma A.l, for to < kl < T, 

WFUxa) C {(p',p) -.peK.XKT, Mp') e K^, ^i^p')) = ^Mp))}- (3.9) 

The definition of - see (3.2) and (3.5) - shows that the set in (3.9) does not intersect 
Diag(T*M). This means that the operator xGr satisfies the hypothesis of Lemma 2.3, 
and by (2.10) we find 

^ I L,^^h-Y,fiE,){xGs-tU,,u,)dtds = Orih^- 

\t-s\>to j 

Combining this with (3.8) and recalling (3.6) and (3.7), we get 

J2 \{Fu„u,)\<{GT-' + OT{h))K 

where C is a constant independent of T and h. By choosing T large and then h small, 
we obtain (3.4). □ 



4. Quantum ergodicity for restrictions 

We will now prove Theorem 1 and we use the notation from the second (semiclassical) 
part of Section 1. To simplify the presentation we put Q = id. The general case is 
similar. 

We start with some geometric observations. The condition (1.5) shows that, in 
the notation of (1.6), Be '■= tteI^e) C T*N, is a smooth manifold with a smooth 
boundary. Any p E Be \ OBe is a regular value of tte] moreover, vr^^(p) = {p+,p_}, 
7Te is a local diffeomorphism near p±, and the involution je is given by 7£;(p±) = p=p. 
The Hamilton vector field Hp is transversal to S^; at p±. 
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To prove Theorem 1 we can assume that [a, b] is a small neighbourhood of a fixed 
energy level E. We then decompose any compactly supported A e "i/l^i^N) as follows: 

J 

A = J2 ^J,e + + (1 - Xe)A, (4.1) 

where 

• Xe e "^^""^^{N) is microlocally equal to id near Be C T*N, and WF/,(Xe) is 
contained in small neighbourhood of Be, 

• Aj^, G ^™°'P(A^), and WFh{Aj^e) is a small open subset of Be \ OBe, 

• G ^™™^(A^), /io- (WF/i(yle)) < e, where /lo- is the symplectic measure. 

The estimate (2.8) in the second part of Lemma 2.2 shows that the contribution of 
(1 — Xe)A is negligible, and that the contribution of A^ will disappear in e — )■ limit. 

Hence we only need to prove Theorem 1 for terms of the form Aj^^. We assume now 
that p G Be \ OBe, and that A G \E'™™''(A^) is microlocalized in a small neighborhood 
V C T*N of p. Choose small 6 > and define the set 

u := {iptix,0 ■■ \t\ < (5,0 e j:E+r nn^l^iv), |r| < 6}. 

If V and 5 are small enough, then we can write U = Ui U U2, where Ui, i = 1,2, are 
open subsets of T*M (one of which is a neighborhood of p+ and the other of p_) and 
moreover, the maps : Ui ^ V x {\t\, |r| < 6}, 

Ki : (pt{x,i)\ — ^ (7rij+^(x,f),t,r) , (x, 1) G n f/^, |r| < 5, (4.2) 

are diffeomorphisms. The maps ki are symplectomorphisms if we consider |r| < 5} 
as a subset of T*]Ri, with r the momentum corresponding to t; in fact, we provide a 
generating function for in (4.6) below. 

Fix a local coordinate system x = {x',Xn) on M such that = {xn = 0}. We 
identify every half-density u{x)\dx\^^'^ on M with the function u{x) and every half- 
density v{x)\dx'\^^'^ on with the function v{x). Consider the operator 

7^ : C°°(M) C°°(iV X Rt), TZu{t) := (e**(^('^)-^)/\)|;v, (4.3) 

then 

hDtTZ = n{P{h) - E). (4.4) 
Take G \&5j°™^(M) microlocalized inside f/^, but such that 

Xi = l microlocally near Kj^(WFh{A) x {|r|, |t| < 5/2})- 
Let x(t) G C^{-5,5) be equal to 1 near [-5/2,5/2]. Then 

5^ := x{t)TZXe : C°°(M) ^ C°°(iV x M^) 
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are compactly microlocalized Fourier integral operators associated to K£. This fol- 
lows from an oscillatory representation of Q''iiPW-E)/h gjyg^ ^if^^ §10.2]. Indeed, in 
coordinates ) and in the notation of [15, Theorem 10.4], 

^^^-^^ ^) ■= f [ e^('^(*'^'°''')-^-^)6,(^, ^, 0, h)u{y) dydr^, (4.5) 



where 

V'(0, x,r]) = X -T], dtijit, X, d^tp) = p{x, d^^)) - E, 

(Pt{x, d^lpit, X, T])) = {dr,^/j{t, X, 7]), T]). 

The microlocalization inside Ue means that d^^p{x,0,C,) 7^ 0, and that implies that 
^'(tx)r{^ is nondegenerate. Hence, 0,77) is a generating function of k^. 

Now, let u be an eigenfunction of P{h) with eigenvalue E' = E + \, where A G 
[-5/2, 5/2]. Then WF?,(u) C p-^E - 5/2, E + 5/2]) and thus 

u\n = (Xi + X2)u\j\f = Biu\t=o + B2u\t=o microlocally near WFh{A). 

Now, by (4.4) each wg := Biu solves hDtWi = Xwi microlocally near WFh{A) for 
\t\ < 5/2. Therefore, we{t) = e'*^/''wi{0) microlocally near WF?,(A) for \t\ < 5/2. 

Take x(^) e €^{-5/2, 5/2) that integrates to 1. Then 

{A{we\t=o),Wk\t=o)L^N) = {(xii) A) we,Wk) l^n xRt) + 0{h'^). 

Therefore, 

2 



{A{u\n), (m|7v))l2(7V) = ^ {A{we\t=o),Wk\t=o)L2{N) + 0{h°°) 

i,k=l 
2 

= Yl ® ^)^^' ^fc)i^(^xMO + C(/i°°) (4.7) 

e,k=i 
2 

= 5^ (Blixit) 8) I)i?^M,M)i2(M) + 0(/i°°). 



We now need to analyse the operators B^i := BKxit) ® A)Bi. This is split into two 
cases. For k = i, B^i is a pseudodifferential operator and Theorem 2 can be applied 
with B = Big. Note that the operator x(t) ^ A is not pseudodifferential, in fact its 
non-semiclassical wavefront set contains points (t, t', x, x') with t = t' and (x, x') in the 
support of the Schwartz kernel of A. However, B^ is pseudodifferential since Bg, B'^ are 
compactly microlocalized and thus we can replace x{t) by a compactly microlocalized 
operator in x{t) ® A, making the latter pseudodifferential. 

We need to compute the symbol of Bu. For that, we use the integral representation 
(4.5) and the stationary phase method, applicable since 8"^^ ^ip is nondegenerate. 
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More precisely, the Schwartz kernel, B^B£{z,y) is given by 

^ ^ >s(V'(i,s,o,r;)-^(t,x,o,c)+.-C-j/-r?)^^^^^ 0, ri)bi{t, X, 0, C) dxdtdridC. 



{2-nh) 



2n 



"^4,4 '^"^'7 ^"^C 

We apply the method of stationary phase in the x,t,7] variables. The stationary point 
is given hj r] = (, {pt{x,0,dx4't(t,x,0,ri)) = {y,ri), and the value of the phase at the 
stationary point is {z — y) ■ The signature of the Hessian is and, by (4.6) and [15, 
Theorem 10.6], the leading part of the symbol in the region of interest is given by 

det dxd^'ip{t, X, 0, () 



\d^Mx,0,dx^{t,S:,0,O)\ 



det d(t^x)dt:^lj{t, x,0,C) 
where t = t{y, () , x = x{y,(), are the critical points. 
Recalling (4.2), it follows that 

a{B;Be) o kJ^ = \d^^p o kJ^ o ttqI'^ near WF^iA) x {|r|, \t\ < 6/2}, 

where ttq : T*{N x R^) T*{N x M^) maps (x, i, t, r) to (x, 0, r). From here and by 
Egorov's Theorem applied to x(t)cg)A, we get a{Bu)oK~[^ = \d^^po kJ^ onQ\^^a{A)xit) 
near {|r| < 6/2}. Then 

/ a{Bu) dfiEi = , "'^i/rnNN / cr(5«) o t^i'^ dxdidt 

J {p=E,} i^j)) J{r=E,^E} 

= , [ i%pr a(i) dx'de, 

where we parametrized T^e^ by (x',^') G Be^ 

Now, we consider the case k ^ £. Then Bke is a Fourier integral operator with the 
canonical transformation Kke := o k£. We want to apply the decorrelation result 
given in Lemma 3.1. 

Using the definition (4.2) of Ke, we see that the canonical transformation n = 
can be described as follows: 

i^{.Vs{.x,l)) = ips{nE'{.x,l)), \s\<6, (x,f) G Sg/ nf/f. 

To apply Lemma 3.1 we need to verify the following: there exists to > such that the 
set 

£:={peU,n ^^t{U,) : 3 t, \t\ > to , vMp)) = ^{Vt{p))} C T*M , (4.9) 

has yUo--measure zero. To see this, suppose that p E £, t is the corresponding time, and 
s,s' G (-6,6) are such that p = ips{x,^), ipt{p) = ips'{x',C), (5,0, (x',f') G J2e' n Ui. 
Then (x',^') = ipt+s-s'{x,^) and the condition ipt{i^{p)) = K^Lptip)) can be rewritten as 

(x,0 G Sij/, (ft+s-s'ix,^) G Sij/, v?t+s-s'(7£;'(5>0) = lE'i^t+s-s'ix,^)). 
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Put to > 26, then t + s — s' 0. It now follows from (1.7) that the set £ from (4.9) 
has measure zero; by Lemma 3.1, the contributions of B^i, k ^ i to the sum (1.9) go 
to as /i —i- 0. 

Going back to (4.7), (A. 4) and (4.8) this means for A satisfying our localization 
assumptions 



Eje[a,b] 



{A^UjIn), (%|iv))L2(JV) - 



I 2 -^^Ej^Ui 



\d^M-'^*EAA)dx'di' 



oil) 



where Vj := fiEj{p~^{Ej)). Assume that the function / used to define restrictions 
of half-densities in (1.8) is equal to a;„, so that u{x)\dx\^^'^ restricts to u{x,0)\dx'\^^'^. 
Using the canonical transformation k,£, we get the symplectic coordinates (x',t,^',r) 
on Ui, in which the measure duEj from (1.10) is equal to V~^\d^^p\^^ dx'd^' . We have 
thus proved (1.9) with / = a;„; the case of general / follows by taking the operator 
|f^2:„/|~^''^^|<9x„/|~^''^ in place of A. This completes the proof of Theorem 1. 



Appendix A. Semiclassical quantum ergodicity with boundaries 

Let (M, g) be a smooth Riemannian manifold with a piecewise smooth boundary, 
dM. That means that M C M where M is manifold without boundary to which g 
extends smoothly, and dM = IJj=i where Nj are smooth embedded hypersurfaces in 
M. Denote by d°M C dM the open set of all points at which the boundary is smooth, 
namely points contained in exactly one of the hypersurfaces A^^ ; the complement dM\ 
d°M has measure zero (with respect to the surface measure on dM). 

We consider an operator P{h) given by (1.1) with Dirichlet boundary conditions. 
(One can take instead any self-adjoint boundary conditions, as long as the Weyl 
law (A. 3) is known to hold.) Let p{x,^) be its principal symbol; we can extend it 
smoothly to T*M. We make the following assumption similar to (1.5), for E e [a, h\. 

X e d°M, V{x) = E ^ dV{x) i NldM. (A.l) 

Then p~^{E) and Tqoj^M intersect transversally. We write 

p-\E) n T^.,,M = Q+UQ-^U Q% 

where (x,^) lies in if the vector HpX G TM is pointing outside of M, in if 
this vector is pointing inside M, and in if HpX is tangent to the boundary of M. 
The covectors in are called glancing, and under the assumption (A.l) this set has 
measure zero inside p-^{E) HT^omM . 

For (x,^) G p^^{E), we define its broken Hamiltonian flow line ipt{x,C,) as follows. 
Assuming without loss of generality that t > 0, we consider the Hamiltonian flow line 
exp(ti7p)(x, ^), defined smoothly on T*M, and let to be the first nonnegative time when 
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exp(tHp){x,^) hits the boundary. If this happens at a non-smooth point of the bound- 
ary (i.e. on dM \ d°M), or if exp(tQHp){x,^) G then the flow cannot be extended 
past t = to- Otherwise, exp{toHp){x,^) G fi^ and there exists unique (xc^o) ^ 
such that the natural projections of exp(tQHp){x,C,) and (xo,^o) onto T*dM are the 
same. We then define (ft inductively, by putting (pt{x,^) = exp(tHp){x,C,) for < t < to 
and (pt{x,$,) = ipt-to{xo,^o) for t > to- For any T > 0, denote by 

}3TCT*Mnp-\[a,b]) (A.2) 

the closed set of all (x,^) such one cannot define the flow (ft{x,$,) on the interval 
[— T, T] using the above procedure. As shown in [14, Lemma 1], for any T the set 
Bt n p^^{E) has measure zero in p~^{E), and for |t| < T, is a volume preserving 
flow on p^^{E) \ Bt- See also [8, p. 310-311] for a symplectically invariant description 
of the broken Hamiltonian flow. Since the flow (ft is well-defined almost everywhere, 
the standard ergodic theory applies to it. 

We will use the following parametrix construction for the Schrodinger propagator 
away from the set Bt- The following lemma is a rephrasing of results of Christianson [2, 
§3.3]: 

Lemma A.l. Fix T > 0. Assume that A G \1'™™^(M°) is supported away from the 
boundary of M and WFh{A) C p-^{[a,b]) \ Bt- Then for each x e C^{M°) and 
for each t G [—T,T], the operator xe~^^^^^A is a Fourier integral operator supported 
away from dM and associated to the restriction of (ft to a neighborhood ofWFh{A) fl 
<y9^^(supp x), plus an 0(/i°^)L2(jv,/)^j;^2(jv/) remainder. The following version of Egorov's 
Theorem holds: 

^^tPI^Ae-^'"-!^ = At,^ + Oih^)L2^M)^LHM), 
where At^^ G \1'^°™p(M°) is supported away from dM and cr(At^^) = x(a o (ft). 

The following basic Weyl law can be proved for the Dirichlet realization of P{h) 
using the standard Dirichlet-Neumann bracketing method (see [9, Chapter 15]): 

(27r/i)"#{j : Ej G [a,b]} = fi^{T*M f] p-\[a,b])) + o{l) as /i ^ 0. (A.3) 

It follows that eigenfunctions cannot on average concentrate near the boundary: 

Lemma A.2. Assume thatx ^ C^(M°) satisfies < x < 1- Then for a' < a < b < b' , 

(2TihY V / (l-x)|Mirc^Vol< / 1 -xt^/Ua + o(l) as /i^O. 

E,e[aM-^^' ir*Mnp-H[a',6']) 

Proof. Take / G C^{a',b') such that < / < 1 and / = 1 near [a,b]. Since / and 
1 — X are nonnegative, it suffices to show that 

(27r/i)"V / il-x)f{Ej)\u,\^dVol= [ {l-x)fip)dfi^ + oil) ash^O. 
■ Jm Jt'-m 
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This holds since the asymptotics for 1 in place of 1 — x follows from (A. 3), while the 
asymptotics for % follows from (2.6). □ 

We can now prove quantum ergodicity for manifolds with boundary: 

Theorem 2. Suppose that (M, g) is a compact manifold with a piecewise smooth bound- 
ary and that uj = Uj{h) are normalized eigenfuctions of the Dirichlet realization of 
P{h). If (1.3) and (A.l) hold, then for any B G '^\{M°) compactly supported away 
from dM , 



0, h^O. (A.4) 



Proof. Take a',b' such that a' < a < b < b' and (1.3) and (A.l) hold for E G [a', 6']. 
(If the flow is no longer ergodic on p^^{E) when E ^ [a, b], we would need to consider 
a', b' close to a, b, for example a' = a — 1/T and b' = b+ 1/T and crudely estimate the 
contribution of [a', b'] \ [a, b] by the Weyl law.) Take large T > and choose a cutoff 
function xt e C^{M°) such that < xt < 1 and 

/ I - Xrdjj,^ <T''^. 

J T* Mnp-^ ([a' ~l,b' +1]) 

Let the function ijj G C^{a' — 1, 6' + 1) satisfy 

V'(-E) / XTdiiE= I cr{B)diJ,E, Ee[a,b']. 

By Lemma A. 2, it is enough to show that for T arbitrarily large but fixed, (A.4) 
holds for the operator B — ilj{P{h))xT, whose symbol integrates to zero on p~^{E) for 
-Eg [a', b']; therefore, without loss of generality we assume that 

/ cx{B) dfiE = 0, E e [a, b']. (A.5) 

By the elliptic estimate (see for instance [5, Proposition 3.2]), we may assume that 
WFh{B) C p-\{a',b')) and in particular B G The set Bt defined in (A.2) is 

closed and has measure zero; therefore, we can write B = B'j, + 5^, where WF/i(5^) n 
Bt = ^ and ||(T(-B^)||i2(p-i[a/j,']) < T^^. By (2.7), the contribution of Bij. to (A.4) goes 
to zero in the limit limr^-oo linisup/j_>0) therefore, we can replace B by B'rp in (A.4). 
Define the quantum averaged operator 



T Jo 
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Then by Lemma A.l, {B'rp)TXT is, up to an 0{h°°) l^^l^ remainder, a pseudodifferential 
operator in v|/™™p compactly supported inside M° and with principal symbol 

t{{B't)tXt) = XtHB'^))t = Y /^^^^^^ ° 
Since each Uj is an eigenvalue of P{h), we can write the left-hand side of (A. 4) as 

Eje[a,b] 

using the Weyl law (A. 3) and Cauchy-Schwarz, we see that it remains to prove that 

limlimsup/i" > || (-By)r'Uj|||2 = 0. 

We can replace (fi'^iT here by {B'j)tXt, as 



limsup(27r/ir || (1 - Xt)^,!!^^ < T"^ 

Ej<^[a,b\ 

by Lemma A. 2. By (2.7), it remains to show that 

lim \\a{{B'j,)TXT)\\L^{p-^[a',b'])) < lim \\{cT{B'j))T\\mp-ma'm = 0- 

1 — ^oo 1 — >oo 

For this, we write for each E G [a', b'], 

\\{a{B'rp))T\\L2{p-HE)) < ll(o-(5))T||L2(p-i(i?)) + II(o-(5t))t||l2(p-1(£;)). 

The first term on the right-hand side converges to when T — )■ oo by (A. 5) and the 
von Neumann ergodic theorem (see for example [15, Theorem 15.1]), while the second 
term is bounded by \\a{B'rj.)\\L2(^p-i(^E))- D 

Appendix B. From semiclassical to high energy asymptotics 

In this appendix we specialize to P{h) = —h^Ag and show how Theorem 1 implies 
the results of [13]. 

Suppose that (M, g) is a compact Riemannian manifold with a piecewise smooth 
boundary in the sense of Appendix A, and with an ergodic broken geodesic flow (ft '■ 
S*M S*M. Suppose that C M is an open smooth hypersurface whose closure is 
disjoint from the boundary. The energy surface p~^{l) = S*M is the cosphere bundle 
of M and Si = S%M is the restriction of S*M to N; Bi = 7ri(Si) = B*N is the coball 
bundle of N and 71 : Si — )■ Si is the reflection across the orthogonal complement of 
the conormal bundle A^*A^ C T^^M. 

The dynamical assumption (1.7) becomes 

The set of p G S^^M satisfying ipt{p) G S'^M 

and y9f(7i(p)) = 'ji((ft(p)) for some t ^ 0, has measure 0. ^ ^ 
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Let {Mj}°^o be the complete set of eigenfunctions of the Laplacian on {M,g): 



-^gUj = XjUj, \\uj\\L2 = 1, = Ao < Ai < A2 < ■ ■ ■ . 

The statement of the theorem uses the standard concept of a (nonsemiclassical) 
pseudodifferential operator on a manifold - see [8, §18.2]. 

Theorem 3. Let N be a smooth open hypersurface satisfying (B.l) with closure disjoint 
from the boundary. Suppose that A G \E'pijg(^) '^^ ^ classical pseudodifferential operator 
on N , compactly supported inside N . Put Vj := uAn- Then 



A" ^ 

A,<A 



0, A ^ cx), (B.2) 



where (j{A) is the principal symbol of A (a homogeneous function of degree on T*N\ 
{0}j, and the measure Vi is defined in (1.11). 

Remark. Theorem 1 allows more general restrictions auj\]s[ + h\~^duUj\Ni for a,b E 
C°^{N). We note also that the nonsemiclassical formulation of quantum ergodicity only 
implies the angular equidistribution of Vj in T*N. That is natural for the standard 
quantum ergodicity since Uj concentrate on S*M but not in this case as f^'s can be 
microsupported anywhere in B*N . That is remedied in the semiclassical Theorem 1. 

Proof. To show how Theorem 3 follows from Theorem 1 we put V = Q and identify 
L^(M, with L^(M, dvolg) by writing half-densities as u{x)\dYo\g |^/^, where u G 
L\M,dMo\g). 

Let X = {x', Xn) be normal geodesic coordinates near (0, 0) G N, in which = 
{xn = 0}, p{x, C,', = + r{x, ^'), and r{x', 0, is the dual of the restriction metric 
g\N- Suppose that / satisfies f\N = 0, \df{x)\g = 1. In the chosen coordinates the last 
condition means that dx^f = 1, f = 0, on N. Hence the restriction of half-densities 
(1.8) obtained using this choice of / shows that, we obtain an identification with the 
restriction of functions u\n G L^(A^, ci vol^i^^). 

We now write out locally the measure z/i from (1.10). In our coordinates, S*M can be 
parametrized by (x', a;„, ^'), r{x',^') < 1, = ±(1 — r{x,^')Y^'^ (the parametrization 
degenerates at r{x, = 1). The Liouville measure is obtained by requiring dfiiAdp = 
dxd^, and 

= = 2(l-r(l,e))V2^^^^' 

In the notation of (1-11) this gives 

where we parametrized S]^M = {{x', 0, ^) : ^1 + r{x', 0, ^') = 1} by (x', G Bi. 
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To pass from the semiclassical result to the special case of the high energy result we 
put Ej = h'^X'j, h = 1/A. The difficulty lies in controlling low frequency contributions 
and estimates (2.8) and (2.9) are crucial for that. 

Let A be a classical pseudodifferential operator of order on A^, with a compactly 
supported Schwartz kernel in N. (Henceforth operators with hats denote polyhomo- 
geneous operators, while operators without hats denote semiclassical operators.) Its 
principal symbol cr{A) is a homogeneous function of degree on T*N. We define 
E ^°(A^) by putting 

A, := Op, {a{A){l - xiWUJ^))), X e C^il 
Theorem 1 shows that for < a < 6, and Vj = Mj Iat, 

■K*a{A,)dvE, 



x{t) 



1, Itl < 1. 



ft" 

h\j(i[a,b] 

We also have 



{A^Vj, t^j)L2(jv,dvol, 



sljv 



0, /i = l/A-^0. (B.3) 



■KVa{A,)dvE, 



7r*a(A) dui + C(e) 



and hence the result will follow once we show that 

J2 \{{A-A,)v„v,)mN4^oi^,^)\ = 0{e) + 0,{h). 

hXj&[a,b] 

Indeed, (B.3) and (B.4) together give, for [a,b] = [1,2], 



(B.4) 



— as /i — )■ 0. 



(B.5) 



{Avj,Vj)L2(N,d^oi^^^)- / 7Tla{A)dui 

hXj€[l,2] J S'^M 

Summing (B.5) for h = 2^A^^, 1 < A; < log2 A, we get (B.2). 
We now prove (B.4). Using (2.9), it will follow from 

limlimsup/i" y \\{A-A,)vj{h)\\^ = 0. 

hXje[a,b] 

For this, we first claim that for any vector field X on N, 

\\iA-A,)hX\\L2^L2<Ce + 0,{h). 

Indeed, the left-hand side of (B.7) is 0{h) if we put an operator in the class ^^phg or 
h"^^^ in place of A — A^, which means that we can reduce to local coordinates, in which 
we can assume X = dy^ and the full symbol of [A — A^)hX in the non-semiclassical 
left quantization becomes, up to ^E'phg + h^h^ terms, 

r{y,m h) := hriia^{y,iri/\7]\){{l-x{\v\)) - (1 - ^(^l^lskA))) • 



(B.6) 



(B.7) 
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However, dyd^r{y,r]; h) = 0{e + h){r])~^'^^ (here the first cutoff gives the 0{h) term, 
while the second cutoff gives the C(e) term); therefore, by the boundedness of 
classical pseudodifferential operators, we get (B.7). 

Now, let Bq G \I'™™^(A^) be a semiclassical pseudodifferential operator equal to the 
identity microlocally near the zero section of T*N, but supported inside an e^/^ sized 
neighborhood of the zero section. Then we can write 

k 

for some vector fields Xk (independent of e) and some B^ e \&^°™''(A^) (with — 
norm (9(e~^/^)); by (B.7), we have 

\\iA-A,){l - Bo)\\l^^l^ < Ct^ + 0,{h) 

and thus by (2.9), the estimate (B.6) holds for [A — A^){1 — Bq). Same estimate holds 
for {A — v4Ji?o, by recalling that \\A — A^\\i2^i2 = 0(1) and using (2.8) together 
with the bound ||(T(i?o)||L2 < 6("~i)/^. This finishes the proof of (B.6) and thus of 
Theorem 3. □ 
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